oo 

o 

O 

(N 



43 

6 



> 
o 

cn 
o 

oo 
o 



X 



Stability of thermal structures with an internal heating source 

Nestor SancheJ*] 
Institute) de Astrofisica de Andalucia, CSIC, Granada, Spain. 

Eugenic- Lopez 
Departamento de Fisica y Matemdticas, 
Universidad Nacional Experimental Francisco de Miranda, Com, Venezuela. 



1 



Abstract 

Abstract: We study the thermal equilibrium and stability of isobaric, spherical structures 
having a radiation source located at its center. The thermal conduction coefficient, external 
heating and cooling rates are represented as power laws of the temperature. The internal heating 
decreases with distance from the source r approximately as ~ exp(— r)/r 2 , being r the optical 
depth. We find that the influence of the radiation source is important only in the central 
region, but its effect is enough to make the system thermally unstable above a certain threshold 
central temperature. This threshold temperature decreases as the internal heating efficiency 
increases, but, otherwise, it does not depend on the structure size. Our results suggest that 
a solar-like star migrating into a diffuse interstellar region may destabilize the surrounding medium. 

Estabilidad de estructuras termicas con una fuente interna de calentamiento. 

Resumen: En este trabajo estudiamos el equilibrio termico y la estabilidad de estructuras 
isobaricas y esfericas que tienen una fuente de radiacion en su centre La conduction termica 
y las tasas de calentamiento externo y enfriamiento son representadas como leyes de potencia de la 
temperatura. El calentamiento interno disminuye con la distancia a la fuente r aproximadamente 
como ~ exp(— r)/r 2 , siendo r la profundidad optica. Los resultados indican que la influencia de la 
fuente de radiacion es importante solamente en la region central, pero su efecto es suficiente para 
que el sistema se vuelva termicamente inestable por encima de una cierta temperatura umbral en 
el centre Esta temperatura umbral disminuye a medida que la eficiencia del calentamiento interno 
aumenta, pero por otro lado la misma no depende del tamaho de la estructura. Nuestros resultados 
sugieren que una estrella de tipo solar que migre a traves del medio interestelar difuso puede llegar 
a desestabilizarlo. 

PACS numbers: 47.50. Gj; 95.30.Lz; 98.38.Am 

Keywords: hydrodynamics; instabilities; interstellar medium 

Descriptores: hidrodinamica; inestabilidades; medio interestelar 
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I. INTRODUCTION 



Most of the structures observed in Astrophysics can be fully explained only if one takes 
into account several physical processes simultaneously, such as dynamics, self-gravitation, 
magnetic fields, radiation transfer, thermal conduction, complex chemical networks, etc. 
Since the seminal paper by Field [ij, a lot of work has been devoted to the study of the 
thermal and mechanical equilibrium of structures under the action of different physical 
mechanisms |2j. One of the motivations behind this kind of studies has been to explain 
the formation of structures via thermal instabilities. Nowadays, it is clearly established 
that thermal instabilities play a major role in several astrophysical contexts, such as solar 
chromosphere and corona, interstellar medium, planetary nebulae, and intergalactic medium 
Q,0, y, 5|, |6j. A full modeling of these systems may be a relatively difficult task requiring 
large computational resources (e.g. see Refs. flflB). On the contrary, the development 
of simple models including few physical processes allows to clarify the specific role of each 
process, setting up the framework for more realistic models. 

When only thermal conduction and a net heat-loss function are considered, a variety of 
isobaric) thermal structures are obtained and stability criteria can be derived analytically 
10]. These criteria can be extended analytically up to the second order [ll|. In spite of the 



simplicity of this approach (many physical processes are not considered) these results can be 
app.ied successful to the atouuc phase of the ISM Q, molecu.ar clouds fl, aud the sola, 
corona [141 ] . Usually, it is assumed that the heating mechanisms are "local" functions, in 
the sense that they depend only on the local thermodynamical state (temperature, density). 
However, two factors may complicate this picture. On the one hand, the heating mechanism 
could be an explicit function of the position in the structure. This situation applies, for 
example, if a heating source is located close to the structure, such as an interstellar region 
with a nearby star (or star cluster), or the solar corona heated by the chromosphere. For 
coronal loops, it has been shown that the existence or not of thermal equilibrium solutions 

n 

and their stability is subordinated to the spatial dependence of the heating [14j . On the other 
hand, if the heating mechanism is associated with radiative processes, the opacity adds a 
non-local character to the heating function making it difficult to perform an analytical study. 



In Ref . 



151 ] the stability of structures heated by radiative processes was analyzed. Generally, 



the results showed that the effect of increasing the opacity is to increase the stability of the 
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thermal structures. However, in order to achieve certain analytical results that work was 
limited to the unrealistic case of a slab-like configuration irradiated simultaneously on both 
sides. 

In the present paper, our attention is focused on investigating the thermal equilibrium 
and the stability of spherical structures heated by a radiation source located at its center. 
Besides the pure academic interest of generalizing the previous results [h], [ll|, 12, [k], Ij]] 
by including an internal heat source, this kind of study could be very useful because this 
is a common configuration, in particular in Astrophysics. We try to keep the approach 
as general as possible and, therefore, our results can in principle be directly applied to any 
situation in which an isobaric, spherical structure has an internal radiation source. However, 
we use values for constants and parameters that are appropriate for typical interstellar 
conditions because in the future we plan to generalize this work in order to be applicable to 
interstellar regions in the Galaxy. It is known, for example, that the diffuse ISM is linearly 
thermally stable but nonlinearly unstable (e.g. .16]). Moreover, it is expected that stars 
systematically pass through the ISM. Thus, it is relevant to study whether these stars can 
sufficiently perturb an interstellar region as to induce thermal instability. Obviously, many 
more physical mechanisms have to be considered to describe this situation adequately but, as 
mentioned before, in a first approximation this study can help to understand the role played 
by this type of realistic heating. In Sec. [Ill we present the basic equations, and in Sec. IHIIwe 
briefly explain the stability criterion we are using. The effect of the internal source on the 
stability is discussed in Sec. [TV] and, finally, the main conclusions are summarized in Sec. [V] 



II. BASIC EQUATIONS 

The equation of energy conservation for a static, thermally conducting fluid, with pressure 
p constant in time t can be written in the form [17j : 

dT 

pcp-fo = v • Hp, t)vt\ + r( P , t) + A( P , t) , (i) 

where T is the temperature, p the density, c p the specific heat per unit mass, and n(p, T) is 
the thermal conduction coefficient. The functions T(p,T) and A(p,T) represent the heating 
and cooling rates per unit volume, respectively. In many astrophysical situations the quanti- 
ties k, T and A can be approximated as being proportional to power laws of the form ~ p a T^ 
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[bsl ]. Under the assumption of isobaricity, we can eliminate the density as an independent 
variable and write 

k = K T k , T = T T m , A = A T n , (2) 

where the constants k , T , A and the indices k, m and n are given by the physical processes 
involved. For instance, a typical interstellar region in the range of temperature 10 2 T < 
10 4 K with thermal conduction by neutral particles, heating proportional to the density, 
and cooled by collisions between particles, may be characterized by k = 1/2, m = — 1 and 
n = -3/2 QQ. 

In this work we consider an additional heating mechanism. Let us assume a spherically 
symmetric region with radius Rf, that has a radiation source located at its center. The 
photons coming from this source can heat the surrounding gas. In the ISM this heating is 
mainly because the UV radiation field is able to photoeject electrons from atoms and dust 



which subsequently thermalize through collisions 
is fl: 



3]. The energy absorbed per unit volume 



r a (r) = ^<jnF(r) , (3) 

where F(r) is the radiation flux coming from the central source, n is the number of atoms 
per cm 3 , a is the photon capture cross section, and the efficiency parameter 7 is the fraction 
of incident energy which is actually absorbed by the gas. The radiation flux decreases with 
distance from the center (r) as 1/r 2 . In order to avoid dealing with divergences at r = 
when integrating Eq. [TJ, we impose a spatial smoothing scale R s . Thus, the radiation flux 
can be represented in the form [20! ]: 

F(0) exp (-r) 
;i + (r/R s 

where r is the total optical depth from the central source to r: 



andr . (5) 

The value of the arbitrary constant a determines how quickly the behavior changes from 
~ F(0) when r <C R s to ~ 1/r 2 when r ^> R s . We have kept a unchanged (a = 10), but 
several tests showed that our results do not depend on a as long as the condition R s <C Rb 
remains fulfilled. The flux at the center can be calculated by assuming emission from a black 
body of effective temperature T : 

F(0) = 7T / B v (T )du , (6) 
Juv 



B v being the Planck's function. The integration is done in the far ultraviolet (900 — 1100 A) 



because this is the range in which photons heat the ISM efficiently 21 
With these considerations, Eq. [1] can be rewritten as 

at r z Or \ or I 



22|]. 



junr yuj t^ y ) / ? \ 



7crnF(0) exp(— t) 
1 + (r/Rs 

This is the integro-differential equation whose stationary solutions (and their stability) we 
want to study in this work. It is convenient, however, to express it in terms of some non- 
dimensional quantities. If T eq is the temperature to which T = A, i.e.: 

T eq = (r /A ) 1/(n - m) , (8) 

then we can define a non-dimensional temperature as: 

9 = T/T eq . (9) 

Given the (external) heating and cooling mechanisms, T eq is a constant representing the 
temperature at the thermal equilibrium when there exists no internal heating (T s = 0). We 
also write the position normalized to the region size: 

z = r/R b ,z s = R s /R b . (10) 

We have kept the smoothing length fixed at the arbitrary value z s = 10 -8 , but, as mentioned 
before, our results do not depend on the exact z s value as long as z s <C 1 is satisfied. The total 
optical depth from the center to the boundary for the particular case T = T eq = constant 
is: 

T o = 7p > I 11 ) 

where the ideal gas equation of state p/ks = nT (being k B the Boltzmann's constant) 
has been used. Therefore, including Eq. [9] to [11] and normalizing the flux at the center to 
the equilibrium flux, i.e. F = F(0)/F eq where F eq is the black body flux calculated using 
T = T eq , we see that the stationary solutions of Eq. [7] are the solutions of: 

1 ,/ ^e k ^-)+x(ff n -0 n ) + 



z 2 dz \ dz / 

rj Fexp [-r / z (l/fl)dz 

d{i + {z/ Zs y*) l/a 



, (12) 



where 

T _of t b 

eq 

and 



KoT k-m+i ' ^ °> 



_ o(p/k B )F eql R 2 h 

To solve Eq. [12] we need two boundary conditions for the temperature. This paper is aimed 
to quantify the effects of an internal heat source on the stability of the steady structures 
taking as reference the previous papers for which the stability criteria are known. Then, for 
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12 



convenience, we assume exactly the same boundary conditions used in Refs. jlO 
Q,Q: 

de , s 

_ = at z = , 15 

dz 

e = e h at z = l . (16) 

The first condition follows from the symmetry of the configuration: the resulting temperature 
distribution from the center z = to the z — +1 side should be identical to the z — — 1 
side. This is the only condition compatible with this property if we impose the additional 
physically motivated requirements that both the temperature and its derivative exist and 
be finite. The second condition, i.e. fixed boundary temperature is necessary in order to 
integrate the energy equation. We write down explicitly this condition (fixed temperature) 
at z — 1 only to be consistent with the previous papers, but actually a fixed value at any 
other position would be equally valid as long as the integration over z can be performed. In 
fact, we are assuming a fixed central temperature and the integration of Eq. [12] is performed 
outwards from this central position (see Sec. II VI) . This is equivalent to converting the two 
point boundary value problem (which requires considerably more numerical effort) into an 
initial value problem. In practice there is no difference between both approaches because the 
result of interst for this work is the entire family of solutions rather than a single solution and, 
therefore, we have to span the whole range of possible central (and boundary) temperatures. 

III. STABILITY CRITERION 

When r] = (no internal heating source) Eq. [12] reduces to the case studied in Ref. 
flol | - In spite of the simplicity of the case r] = 0, its study has already led to a better 
understanding of the general problem of stability of structures, even allowing to extend the 
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FIG. 1: The central temperature (Oq) as a function of the boundary temperature (9b) with k = 1/2, 
m = —1, n = —3/2, for the case r\ = and the labeled values of A. 



study analytically up to the second order 



We see that for the case rj = the trivial 



solution 6 = 1, i.e. the thermal equilibrium solution T = T eq , is solution o 



been shown that this solution is linearly unstable if A fulfills the condition lo| : 

_2 



Eq. 



A > A r 



7T 



It has 



(17) 



m — n 

The dimensionless parameter A (Eq. [T3l) measures the ratio between heating (or cooling) at 
equilibrium and thermal difussion through the scale length R^. Moreover, thermal diffusion 
is a stabilizing factor because it tends to diminish any temperature fluctuation. Thus, 
physically, the condition A > A cr j means that the generation of heat is so large that it cannot 
be removed efficiently by thermal difussion, and therefore a thermal instability develops. The 
analysis of the rest of the stationary solutions (different from the trivial one) cannot be done 
analytically but the stability can be inferred from the position on the curve 6 vs. 6\>. In 
Fig. [T]we show the plot resulting from integrating Eq. [12] with 77 = 0. The trivial solution 
#0 = 0b = 1 (shown as a filled circle) is located on the positive slope branch of the #0 vs. 0& 
curve for the case A < X cri . This case corresponds to the stable solution according to Eq. [13 
On the contrary, the trivial solution is on the negative slope branch for the case A > A cr j 



S 



(unstable solution). Moreover, for A = A cr j (marginally stable solution) the trivial solution 
separates the positive and negative slope branches. According to Ref. n]] (see also Refs. 
17, Q, S]); the positive slope branch can be identified as the stable branch. Therefore, 
stationary solutions located on the positive (negative) slopes branches are expected to be 
linearly stable (unstable). This is the stability criterion we are using in this work. 

For the more general case (77 7^ 0) we see that 9 = 1 is not a stationary solution of 
Eq. [121 Both the spatial dependence of the internal heating and its non-local character 
(i.e., the local heating depends on material state in other points through the optical depth) 
complicate the mathematical treatment in such a way that it prevents us from finding an 
analytical stationary solution. Thus, a stability analysis is not possible but nevertheless we 
expect the above stability criterion to remain valid. The argument for this claim comes from 
the relative behaviors of thermal conduction and net heating/cooling with respect to the 
temperature. When two stationary solutions are obtained for a given boundary temperature 
(for example when 6^ > 1 for the case A = A cr j in Figured]), only one of them is stable (see 



Refs. 



17, 



231] ). Let us assume, for example, that both thermal conduction Q(T) and net 



heating rate L(T) are increasing functions of temperature. Two stationary solutions means 
that Q = L at two temperature values (T = T eq ). Necessarily, we have dQ/dT > dL/dT 
at one point and the opposite in the other. The point at which dQ/dT > dL/dT is stable 
because any small temperature variation will be eventually damped by thermal conduction. 
On the other hand, if the heating rate increases with temperature faster than the thermal 
conduction then any small temperature rise will grow. This argument should remain true as 
long as the temperature dependences remain the same. For a fixed position r, the dependence 
of the internal heating on T is the same as in the previous papers (~ T , i.e. proportional 
to density). This is the reason why we expect the stability criterion to be valid for the case 
under study. 

In any case, we performed some numerical tests in order to check the validity of this 
criterion for the more general and complicated case of a non-local heating mechanism. For a 
given set of free parameters (A and rj) we first find numerically the corresponding stationary 
solutions for central temperatures around the "turning point temperature", i.e. around the 
central temperature value that separates the positive and negative slope branches. Then we 
impose different temperature perturbations (we use gaussian, exponential, and sinusoidal 
profiles) around randomly chosen positions in the structure. What we did was to check 
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that dT/dt changes its sign exactly at the turning point, in such a way that on the positive 
slope branch the sign is opposite to the perturbation (stable branch) and viceversa. We 
performed many tests with different combinations of parameters and the results obtained 
were the same. Notwithstanding this is not a demonstration, these tests allow us to be much 
more confident about the validity of the stability criterion for the case under consideration. 



IV. EFFECT OF THE INTERNAL HEATING 



To numerically solve Eq.[T2]with the boundary conditions (11511161) we use the Runge-Kutta 



method of order 4 with adaptive step size control [24| . The distinctive feature of Eq. [12] is the 
fact that to evaluate the internal heating term we need to know the solution we are looking 
for. The density is linked to the temperature (because of the isobaric assumption) and 
therefore the integral in Eq. ultimately depends on the temperature distribution from the 
source (as we have indicated explicitly in Eq. fT2i) . Thus, we have to proceed by successively 
iterating Eq. [12] until 9{r) converges to the stationary distribution. Then, given a central 
temperature 8q, i n the first step we assume an initial distribution 9{z) = 8 to obtain an 
initial estimate of the optical depth r(z). Then we use this result to numerically integrate 
Eq. [12] and obtain a new temperature profile. Now this temperature distribution is used 
to recalculate t(z), and so on. The initial conditions are always satisfied, thus the central 
temperature remains constant during these steps. The procedure is then repeated until a 
convergence criterion is satisfied. In particular, we require that temperature values computed 
at two successive steps differ less than a prefixed small tolerance. Furthermore, we require 
this condition to be met at several points inside the structure. The final result is the whole 
temperature distribution from the center to the boundary. The central temperature can then 
be changed to search for other solutions with their corresponding boundary temperatures. 

The coefficients A and i] are determined by the physical processes involved (pressure, 
external heating, cooling, thermal conduction, chemical composition, etc.). Here we consider 
-ypical ISM conditions. The mean pressure of the warm gas in the ISM isp/ks — 10 3 K cm~ 3 



251 ] , which is heated mainly (in absence of nearby stars) by the galactic radiation field. Under 



these conditions, we have Tq = 2.7 x 10 23 , Aq = 10 21 , m = —1, and n = —1.5 3], [26] 



m 



this work T and A are in units of erg cm 3 s 1 ). This yields an equilibrium temperature 
T eq ~ 1370 K. Under these conditions, the thermal conduction is mainly due to neutral 
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FIG. 2: The distributions of temperature with k = 1/2, m = —1, n = —3/2, for the case Rb = R cr i 
and 7 = 10~ 3 . Dashed line corresponds to the highest central temperature shown (8$ = 4.4). 

particles and therefore Kq = 2.5 x 10 3 and k = 1/2 j^]. The cross section is assumed to be 
a = 6.3 x 10~ 18 cm 2 |3j. We keep as constants these fiducial values whereas the region size 
(Rb) and the internal heating efficiency (7) are treated as free parameters. From Eqs. [13] 
and [T7] we get that the case in which the trivial solution is marginally stable corresponds 

1/2 

to structures with radius R rr i = (\ cri k T e fc - m+1 /r ) ~ 0.12 pc. The total optical depth 
(in thermal equilibrium) would be tq ~ 1.7. As mentioned before, A measures the ratio 
between net external heating at equilibrium and thermal conduction through the structure. 
The trivial solutions is stable when A < A cr j because thermal conduction smooths the energy 
gradient more efficiently in relatively small structures (Rj, < R cr i)- On the other hand, rj 
gives in some way the relative importance of the internal heating and thermal conduction 
terms. Thus, given three different mechanisms of energy transport (thermal conduction, net 
external heating, and internal heating) there are only two free parameters (Rb and 7, related 
to A and 77 through Eqs. [T3land [T4l determining the relative importance among them. 

The temperature distributions inside the structure for the case Rb = R cr i and 7 = 10~ 3 
[19I are shown in Fig. [2j Each line in this figure is one solution 6(z) of Eq. [12] for a given 
value of central temperature 8q and, therefore, of internal heating. We see that there are no 
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solutions with boundary temperature (6b) below 1. Then, all the solutions with #o < 1 have 
positive temperature gradients. As we will see later, these solutions are thermally unstable. 
For 6*o > 1, we have that 6 increses as 0& increases, but for 6 > 4 this behavior changes 
notoriously. For the highest central temperature shown, we can see a sign change in the 
temperature gradient around z ~ 0.55. To understand this behavior, we have calculated the 
relative contribution of each term in Eq. [12] and we have seen that at central temperatures 
#o~2 the internal heating is always much smaller than the (external)heating/cooling term. 
Thus, the situation is almost the same as if 7 = 0, i.e., (external)heating/cooling balanced 
only by thermal conduction. In contrast, at higher #0 values the internal heating dominates 
over the (external)heating/cooling, but this internal heating decreases very quickly with the 
distance from the center and then a change in regime occurs at some point (depending on the 
central temperature). For the highest central temperature shown in Fig. [2l this change in 
regime occurs precisely at z ~ 0.55. We will see in next figure that this central temperature 
is very close to the maximum central temperature for which there exists a solution. Above 
this maximum, the temperature gradient becomes so large that the temperature goes down 
to unphysical negative values (i.e. there is no physical solution). For central temperatures 
between ~ 2 and this maximum temperature the internal heating term always dominates at 
z ^ 0.5. Therefore, the influence of the radiation source seems to be important only around 
the central region (for the values of constants and parameters considered here) , but its effect 
can be enough to make the system unstable, as we will show next. 

The 60 vs. 6b diagram for Rb = R cr i and various values of 7 is shown in Fig. [3j Each 
point in this figure represents one solution of Eq. Q2] (i.e., one line in Fig. [2]), and each 
line represents a family of possible solutions for each pair of values (Rb, 7). Note that all 
curves cross through the trivial solution (6q = 6b = 1, shown as a filled circle). This is 
contrary to the expected behavior because 6 = 1 = constant is not a solution of the energy 
equation. But, as mentioned before, at central temperatures around 6q = 1 the internal 
heating can be neglected and the situation is equivalent to the case 7 = 0. According to the 
stability criterion (Sec. IIIII) . the solutions located on the positive slope branch are stable, 
and viceversa. For the case 7 = we recover the solutions studied in Ref. n]], where the 
trivial solution is marginally stable and the other solutions with 6 > 1 are stable. For the 
cases 7 ^ Owe observe that there exists a maximum central temperature (hereinafter 6o^ max ) 
above which the solutions become thermally unstable. For the case 7 = 10 -3 this threshold 
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FIG. 3: The central temperature (9q) a s a function of the boundary temperature (9b) with k = 
—1/2, m = —1, ra = —3/2, for the case Rb = R cr i and the labeled values of 7. The thermal 
equilibrium solution is indicated by a filled circle. 

occurs at 9o, max ~ 3.9, but clearly 9o, ma x decreases as 7 increases. 

In Fig. S] we have fixed the internal heating efficiency (7 = 10~ 3 ) and we have varied 
the radius Rb- At low central temperatures the behavior is similar to the case 7 = 0, i.e., 
the trivial solution is stable for Rb < R cr i and unstable for Rb > R c h- At relatively high 
central temperatures, we obtain that the threshold temperature 6 0tTnax for thermal instability 
(indicated as open circles in Fig. HJ) is approximately constant. This interesting result is a 
direct consequence of the "local" character of the internal heating: this heating mechanism 
decreases too quickly with z that its destabilizing effect does not depend on the structure size. 
It can be seen that there is a small variation in 6o ymax for the extreme case Rb = 10~ 6 R cr i, 
but this value is so small (~ 5-R sun for our fiducial values) that its validity is questionable. 

In the case of a star that migrates into a diffuse interstellar region, the quantity 0max 
would represent the photospheric (effective) temperature of a star above which the surround- 
ing medium becomes thermally unstable. We have to point out that this simple interpreta- 
tion is possible only if the time-scale for thermal conduction is shorter than any other time 
scale (e.g. the dynamical time) that could be involved in the problem. For our fiducial case 
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FIG. 4: As in Fig.[3]but for the case 7 = 10 3 and the labeled values of Rb- The thermal equilibrium 
solution is marked by a filled circle, and the open circles indicate the threshold solution 8 = #o,maz • 

with 7 = 10~ 3 and Rb = R cr i — 0.12 pc we obtain 8 0tmax ~ 3.9. An interesting result is that 
8o,max does not depend on, or depends very weakly on, Rb (at least under the conditions and 
assumptions considered in this work). In other words, our results suggest that a solar-like 
star with surface temperature ~ 5300 K may possibly destabilize the medium in its vecinity, 
almost independently of the region size. After this, the surrounding medium will undergo 
a transition to another state (or phase). Differences in the chemical composition and/or 
the amount of dust from region to region in the Galaxy (or in other galaxies) may change 



the value of 7 19j, and consequently the maximum central temperature. Here we obtain 
0o,max = 3.6 for 7 = 10~ 2 and 9 , ma x = 3.2 for 7 = 10" 1 . 

If stars eventually migrate into interstellar regions, then thermal instabilities caused by 
these stars may play a non-negligible role in stimulating phase transitions in the ISM. How- 
ever, the present model is clearly too simple to be directly applicable to the ISM. The 
strongest restriction is perhaps isobaricity, which is a difficult condition to be fulfilled even 



in the diffuse ISM 



251 ]. Notwithstanding these limitations, our results suggest that thermal 



instabilities driven by stars could be a key mechanism in the ISM. 
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V. CONCLUSIONS 



We have studied the equilibrium and stability of isobaric, spherical structures with an 
internal radiation source located at its center. We find that there exists a threshold central 
temperature above which the structure becomes thermally unstable, and this temperature 
decreases as the internal heating efficiency increases. Interestingly, the threshold temper- 
ature does not depend on the region size, because the destabilizing mechanism has only a 
local effect. We have shown that thermal instability triggered by an embedded source may 
have interesting consequences in the ISM evolution which should be addressed in future, 
more complete, studies. 
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